#!/usr/bin/env python3
"""Neuroinflammation Analysis Script
Author: Juan Moises de la Serna
"""
import pandas as pd
import seaborn as sns
import matplotlib.pyplot as plt

df = pd.read_csv('biomarkers_extended.csv')
print(df.groupby('diagnosis')[['csf_abeta42_mean','csf_tau_mean','plasma_nfl_mean']].mean())

sns.boxplot(data=df, x='diagnosis', y='plasma_nfl_mean')
plt.title('Plasma NfL by Diagnosis')
plt.tight_layout()
plt.savefig('nfl_by_diagnosis.png', dpi=150)
